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A matterless double slit 

Ben King^, Antonino Di Piazza^ and Christoph H. Keitel^ 
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Double-slits provide incoming photons with a choice. Those that survive the passage have 



^ ■ chosen from two possible paths which interfere to distribute them in a wave-like manner. 
§ ! Such wave-particle dualitj^' continues to be challenged^ and investigated in a broad range of 
^ ; disciplines with electrons^, neutrons', helium atom^, Ceo fuUerenes', Bose-Einstein condensa- 
tei^and biological molecules^^ All variants have hitherto involved material constituents. We 
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p |! present a matterless double-slit scenario in which photons generated from virtual electron- 

q 

^ ■ positron pair annihilation in head-on collisions of a probe laser field with two ultra-intense 

>>. laser beams form a double-slit interference pattern. Such electromagnetic fields are pre- 

dieted to induce material-like behaviour in the vacuum, supporting elastic scattering be- 

> ■ tween photom^EHSl. Our double-slit scenario presents on the one hand a realisable method 
oo 
cn 

O ! to observe photon-photon scattering, and demonstrates on the other, the possibility of both 



controlling light with light and non-locally investigating features of the quantum vacuum's 
structure. 



?H ■ According to both special relativity and Heisenberg's uncertainty principle, virtual electron- 

positron pairs spontaneously pop into and out of existence in vacuum, on a time scale too short 
to leave a trace. However, it is the polarisation of these pairs under an applied electromagnetic 
field which is predicted to provide a rich variety of non-linear processed. A fundamental scale 
for such vacuum polarisation effects is set by the critical field of quantum electrodynamics Ecr = 
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\f\^w?c' jfie = 1.3 X 10^^ Vcm~^, for electron mass m and absolute charge e (in our units the 
fine- structure constant reads a = /A-kHc ~ 1/137), corresponding to a laser intensity of Icr = 
2.3 X 10^^ Wcm"^. An electric field of this order is strong enough to provide a virtual electron- 
positron pair an energy equal to its rest energy 2mc^ in the fleetingly short time h/mc^ ~ 10"^^ s 
in which the virtual pair "lives," promoting it to reality before the individual particles eventually 
annihilate with one another. Even at much lower intensities / such as provided by "strong" or ultra- 
intense (/ > 10^^ Wcm^^) laser fields, the polarised vacuum is predicted to exhibit birefringence 
and dichroism'^, to cause photons to "merge" or to "split" and even allow them to scatter, all 
of which awaits experimental confirmation in laser fields. Recent advancements and proposals 
for the upcoming ELl^^ and HIPEFP laser facilities demonstrate a maturing of a technology that 
will supply intensities of the order 10^^-10^^ Wcm^^, which are sufficiently high to test quantum 
electrodynamics in this relatively unprobed regime. 

When driven by a strong electromagnetic field, the virtual electron-positron pairs generate a 
polarisation and a magnetisation in the vacuunJ^^ (see Supplementary Information): 

P(t,r) = —-J2{E'-B')E + 7{E.B)B] (1) 
M(t,r) = -^^[2{E'-B')B-7{E-B)E], (2) 

for electric and magnetic fields E(t, r), B(t, r). From these expressions, we can form the useful 
analogy of the polarised vacuum as a solid with non-linear response, which, instead of comprising 
tangible dipoles, hosts transient polarised virtual parti cle-antiparticle pairs of dimension approxi- 
mately equal to the Compton wavelength Ac = h/mc ~ 10^^^ cm. These evanescent pairs mediate 
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a non-linear interaction between fields which becomes more significant the larger the ratio of the 
applied to the critical field becomes. When discussing strong electromagnetic fields, we are thus 
referring to a regime totally forbidden in classical physics in which the linear superposition prin- 
ciple in vacuum no more applies. 

Taking the solid-state paradigm one step further, using an ultra-intense laser split into two 
beams, the vacuum can be "activated" by polarising two slit- like regions (see Fig. 1). When these 
regions are probed with a second, (almost) counter-propagating laser, one can imagine creating a 
real photon-photon double-slit experiment. This employs Babinet's principle, which states that the 
diffraction pattern of an aperture is identical to that of an opaque obstacle with the same shape as 
the aperture, justifying our labelling of the two polarised regions as "slits," although they are actu- 
ally the material-like portion of the scenario. Since accelerated charges radiate, when the polarised 
vacuum is agitated by the applied field, it forms a source or vacuum current of electromagnetic 
waves, Jvac(i!:,r). The modified wave equation incorporating vacuum polarisation effects reads 
(see Supplementary Information): 

V'E-i^E = J™,(t,r), (3) 

where Jvac(t, r) = (l/c)V A {dtM) + {l/c^)d'^V - V(V ■ P). This current is then responsible 
for the generation of two fields ,(/!:, r) with i = 1,2 each in the centre of the two slits. These 
fields then interfere to produce the characteristic double-slit diffraction pattern. In Young's origi- 
nal experiment all other incident light was stopped, whereas in our scenario, the probe laser can 
form a dominant background. Exploiting the wide extension of the field generated in the slits in 
comparison to the relatively tight focusing of the probe field, allows us to consider regions on the 

3 



detector plane where the probe background is effectively negligible. This is a significant advance- 
ment with respect to previous calculationi'^, which although based upon the same fundamental 
physics of quantum electrodynamics and quantum vacuum fluctuations, focused only on the ellip- 
ticity and rotation of the polarisation direction acquired by an X-ray probe when it collides with 
a single strong optical standing wave. By introducing the double-slit, the interference pattern of 
photons generated in the annihilation of virtual electron-positron pairs occurring at different points 
in space then becomes a useful, measurable quantity and provides in principle both non-local in- 
formation about the vacuum current and insight on the wave-particle duality of vacuum- generated 
photons. 



With regard to the corresponding experimental implementation, it is pertinent to consider the time- 
averaged total signal /^(r) on a detector plate whose origin is situated in the far field at r = (0, y, 0), 
comprising the scattered fields Ed(t, r) = E^^ r) + ^d,2it, r) and the unperturbed probe field 
Ep(t,r): /i(r) = (|Ep(t, r) + E,(t, r)p) = /^(r) + /p,(r) + /,(r). Here /^(r) = (Ep(t,r)2), 
Ipd{r) = 2(Ep(t, r) ■ Ed(t, r)) and /^(r) = {Ed{t, rf), with () denoting a time average. 

In terms of apparatus, the probe laser should be optical in order that the diffraction pattern 
is sufficiently large and resolvable. We consider the following nowadays easily-obtainable param- 
eters of 100 fs pulse duration, intensity 4 x lO^^Wcm^^ and wavelength Ap = 527 nm (achiev- 
able using the second harmonic of readily-available 1054 nm lasers with an intensity attenuation 
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Figure 1: A matterless double-slit set-up. Two ultra- intense Gaussian laser pulses with wavevec- 
tors ks 1 and 2 are tightly focused by the lenses Li and L2 (almost) antiparallel to a probe beam 
with wavevector kp of much wider spot radius (see also inset a)). The vacuum current, activated in 
the interaction regions of the probe and the strong laser fields, generates photons which interfere 
to produce a diffraction pattern on the screen S. The screen is placed between the focusing mirrors 
at a distance y along the propagation axis of the probe from the interaction centre and has a hole 
in the centre allowing the probe beam to pass undisturbed (see also inset b)). The directions of the 
spatial co-ordinates x,y,z and the angle 9 between the strong field and the probe field Ep are 
defined in the inset a). 
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of around 2.6). For the strong-field laser we anticipate parameters available from the upcoming 
ELI and HiPER facilities of intensity 5 x lO^'^Wcm^^, pulse duration = 30 fs, wavelength 
Xs = 0.8 nm, spot radius w^^o = Xg = 0.8 fim (corresponding to a laser power of 100 PW) and 
laser repetition rate 1 min. 

The success of the solid-state perspective is then displayed by plotting the bare vacuum 
signal /d(r) and observing the accuracy of the famous double-slit formula {n + 1/2) Ap = Dsini!) 
for predicting minima, indicated by crosses on the x-axis in Fig. 2a. D is the distance between 
the centres of the two ultra-intense lasers, ?9 = tan"^(x/y) for detector distance y along the axis 
from the interaction centre and transverse displacement of the minimum on the detector x, with 
different integers n corresponding to different minima positions. We have chosen to separate the 
strong beams by -D = 80 Ws,o = 64 fim and to polarise the probe at = n/2, focused onto a spot 
of radius Wp o = 290 yum and to place the detector at y = 5 m. The choice of D is sufficiently 
large such that the diffraction pattern can be observable also with more realistic (broader) strong 
beam transverse intensity profiles The single-slit limit of zero strong-field beam separation is 
depicted in Fig. 2b, in which all fringes have disappeared. The corresponding diffraction pattern 
does not show typical diffraction rings due to the "slit" having edges that are not sharp. 



Turning to quantitative results, we envisage verifying the phenomenon by either a full or partial res- 
olution of the interference pattern, or else by simple counting of diffracted photons. At our probe's 
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Figure 2: Light-light diffraction pattern. For the experimental parameters given in the text, the 

logarithmic plot of the vacuum signal in a) reveals a series of bright and dark fringes resembling the 

characteristic double-slit pattern. The prediction of the classic {n + l/2)Ap = D sin ^9 formula for 

minima is indicated with crosses plotted on the x-axis. For a radius greater than 1.5 cm, indicated 

on the figure as a black contour, /^(r) is much greater than the /p(r) and /pd(r) background. For 

the same parameters as a) but at zero strong-beam separation, the logarithmic plot of the vacuum 
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signal in the far- field detector plane is shown in b). 



wavelength of 527 nm back-illuminated CCDs (charge-coupled devices) have an efficiency of 90% 
From numerical results for the aforementioned typical experimental parameters, on a region in 
which /d(r) is more than one-hundred times larger than Ip{r) and Ipd{r), taking into account CCD 
efficiency, we expect per shot of the strong field, 4 photons from the vacuum signal. By repeating 
the experiment first in the absence of the strong beams and then the probe and vice versa, one is in 
principle able to account for possible background photons coming from those beams. Background 
photons with a frequency different to that of the probe could be excluded by placing frequency 
filters in front of the detector screen. The presence of a thermal photon background can then be 
neglected operating at a typical temperature of the order of 300 K. Moreover, a good vacuum 
quality of the order of 10^^-10^^ torr is required in order to neglect diffraction effects due to the 
presence of residual gas in the interaction region. We have also ensured that with the above numer- 
ical parameters, alterations to the vacuum signal due to the pulse shape of the strong beams can be 
consistently neglected. 

One can form the visibility of the diffraction pattern (see 1^ and Supplementary Information), 
to determine how many photons are required before fringes can be adequately differentiated. For 
a scenario in which a 15 cm x 15 cm CCD with a central circular aperture of 1.5 cm radius, placed 
as indicated in Fig. 1, detects vacuum signal photons for the aforementioned experimental param- 
eters, a theoretical maximum visibility of 47.6% can be reached. After modelling experimental 
trials numerically, it was found that ~ 1000 photons were required before the statistical fluctua- 
tions around this analytical value were reduced to less than 10%, corresponding to an operating 
time of approximately four hours (see Supplementary Information). 
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By exploiting the polarised vacuum in such a scenario with lasers available in the next few 
years, one can take Young's famous experiment one step further and create a truly quantum double- 
slit set-up comprising entirely of light. In addition, by counting photons or measuring the intensity 
pattern directly, such a method can be employed to probe the quantum vacuum and to study its 
structure as predicted by quantum electrodynamics. 
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Supplementary Information 



In the limit of electromagnetic fields E(t, r) and B(t, r) with amplitude much less than the critical 
fields E„ = VA^m^cyhe = 1.3 x lO^^ V/cm and 5„ = Vl^m'^cyhe = 4.4 x 10^^ G (in our 
units the fine-structure contant reads a = e^/Anhc ~ 1/137) and with wavelength much larger 
than the Compton wavelength Ac = h/mc = 3.9 x 10^^^ cm the vacuum Lagrangian density of 
the electromagnetic field including quantum correcting terms due to vacuum polarization is given 

^=liE'-B') + ^, [{E' - Br + 7(E • B)^] , (4) 

where units h = c = 1 axe employed and terms proportional to represent quantum correc- 
tions much smaller than the Maxwell Lagrangian density (E"^ — B'^)/2. In our scenario, the total 
electromagnetic field consists of two strong Gaussian-focused beams linearly polarized along the 
x-direction that propagate along the negative y-direction antiparallel to a weaker Gaussian-focused 
probe beam linearly polarized at an angle 6 to the x-axis. The strong beams have frequency uj^ 
(wavelength As = 27r/a;s), peak electric field and are centred at (x, z) = ±(xo, ^o) with 

waists w^{y) = w^^^Jl + {y/y^^^Y and Rayleigh length ?/r,s = ^sw'^fl/'^- The probe beam has 
frequency tUp (wavelength Ap = 27r/wp), peak electric field Ep^ and is centred at (x, z) = (0, 0) 
with waist Wp{y) = Wpfl^/l + {y/y^^pY, Wp^ > w^^ and Rayleigh length y,^p = cUpWp o/2. Here 
we concentrate on the currently unknown leading order contribution of elastic real photon-photon 
scattering and in our analysis correspondingly take a solution of Maxwell's equations to first order 
in (w^j.o/z/rj) with j = s, p (see e. g. I^for more details on the approximation used here). We also 
neglect the angle between the strong laser beams (see Fig. 1 in the main text) and assume they 
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propagate in the same direction. Angles of the order of 0.1-0.2 rad can be in principle achieved ex- 
perimentally and, following our numerical simulations, lead to corrections of the order of 10-20 %, 
respectively. The following equations are used to represent the system: 



E(t,r 
Es(t,r 
Es,i/2{t,r 
Ep(t,r 



where 



and 



E,(t,r) + Ep(t,r) 
[Es,i(t,r) + Es,2(t,r)]x 

-^£s,o{x =F xo, y,zT zq) sin (uj,{t + y) - /s(x =F xq, y,zT zq] 

Ep{t, r) [x cos 6* + z sin 9] 

Spfi{x, y, z) sin {ujp{t - y) + fp{x, y, z)^ , 



(5) 
(6) 



(7) 



^j,o{x,y,z) 



-{x^+z^)/w]{y) 



(8) 



fj{x, y, z) = tpj + tan ^ (—) 



Ujy x"^ + z"^ 



(9) 



.y,j/ 2 y^ + ylj 

where is a constant phase offset, j = s, p. The inclusion here of defocusing terms in the probe 
field Ep(t, r) scaling as y/?/r,p is a significant improvement on previous result^, allowing us to 
investigate the vacuum polarization effects also in the so-called far region, where the observation 
distance y is so large that y/yr,p S> 1. This is essential here, as the suggested experimental setup 
requires the observation screen to be located far from the interaction region (in the numerical 
example considered in the main text we have y/yr,p ~ 10). 



By applying the principle of least action to the Lagrangian density in Eq. dH), one obtains the 
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wave equation for the total electric field: 

V'E-S^E = J^c, (10) 
where the "vacuum current" Jvac(^5 r) can be written as 

Jvac = VAdtM- V(V ■ P) + a^P, (11) 

with P(t, r) and M(t, r) being the vacuum polarization and magnetization respectively: 

P := ^[2(i?2-i?2)B + 7(E-B)B] (12) 
M := -^[2(E^-i?^)B-7(E.B)E]. (13) 

The wave equation (fTOl) can be written formally as the integral equation 

E(t, r) = Eei(t, r) + y dt'dh'D{t -t',r- r') Jvac(t', r'), (14) 

by employing the Green's function D{t, r) = —1/(27?)"^ J dujdkex-p[—i(ujt — k ■ r)]/ (w^ — |kp) 
(see, for example,!^. The first term in this equation is the classical solution that in our case is given 
by Eci(t, r) = Es(t, r) + Ep(t, r). The second term arises due to the quantum interaction between 
the probe and the strong fields, which we label the diffracted field Ed(t, r), and is calculated by 
substituting the zero-order solution Eci(t, r) into the vacuum current Jvac(^5 r). Since our probe 
and strong fields are monochromatic it is convenient to work in the Fourier-transformed frequency 
space and the diffracted field can be written as: 

gj(a;p(t-r)+i/ip) 

Eait, r) = Ed(r) + c.c. . (15) 

The Fourier amplitude Ed(r) is then given by 

Mr) = 'f'^tMr), (16) 

fc=l 
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where 1^^ = is the strong field intensity, v is the polarization vector with: 

^ 4(l + f)cos^ ^ 

-(1 - f)(4f cos^ + 7f sin^) 
7(1 + f) sine 
and Xk(r) are the four integration volumes: 



V 



o 



X\2 /;ZN2 

r J \r 



(17) 



3^/ 



d r 



'l + {y'/yr,.y),/T+WKv 



(18) 



, x'^ + y'^ + z'"^ xx' + yy' + zz' 



(xx' + y?/' + zz'Y 
2r 

+ i tan" 



-(x +z +Xo + 2:o) 



x'2 + z''^ 



^ y' ioopy' x'^ + z'"^ 



wl{y') 



2 y'^ + yl 



+ (x'xo + z'zq) 



4/3k , 2iT^uj,y' 



+ 



(19) 



where /3i = l,/32 = -1 and (3^ = 13^ = 0, Ti = Ta = 0, = 1 and 14 = -1, where 
A^/^s = z^s,2 — V's.i- I V ■ v| = v^ is maximized for cos 26 = —1, and without loss of generality, we 
set A^s and zq to zero. 

Since the diffracted field contains spatial integrals over the probe and the strong fields, its 
decay length in the transverse x-z plane results in being much larger than that of the probe field. 
It can be shown that the integrals I^ir) and X4(r) are negligible with respect to Xi(r) and X2(r), 
which depend only on the physical parameters of one of the strong beams respectively and therefore 
describe the interaction of the probe field with each "slit." Therefore the diffracted field amplitude 
E(i(r) can be written as Ed(r) = Ed,i(r)+Ed,2(r) with the subscripts 1, 2 referring to the respective 
terms in Eq. (fT6] ) and the quantities Ed,t(t, r) (i = 1, 2) employed in the text derived from Eq. (fTST ). 
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The analytical expression for /d(r) = (Ed(t, r)^), with () denoting a time average, in the 
limit of no probe focusing, y^^p — )• oo, Wpfi ^ Ws,o and x/r,z/r ^ y/r ^ 1 and zero beam 
separation xq = zo = was also derived and found to have excellent agreement with numerical 
results. As a further check, we derived the ellipticity e induced in the probe and compare this to 
the expression for two, parallel-propagating, colliding lasers in the refractive-index (i. e. short 
observation distances, y 0), crossed-field (Ws 0) limit found in other literature^^. One 
achieves the result e{6 = 7r/4) = (2a;7r/15)(/s,o/-^cr)(^y/Ap), where ly is the effective interaction 
length of the two sets of beams, ly = 7r|/r,p2/r,s/ (2/r,p + l/r,s), which then agrees in the limit y^^p — )■ oo 
with the literature!^. 

For a fixed y, we can approximately maximize the region in which /d(r) ^ -^p(r) + -^pd(r)> 
with /p(r) = (Ep(t,r)^) and /pd(r) = 2(Ep(t,r) ■ Ed(t, r)), and hence maximize the signal-to- 
noise ratio by ensuring the fastest decay of the probe and mix-term background in this far-field 
plane. This occurs when the Gaussian variance is minimized, corresponding to a probe focused to 
a waist of Wp ^ = ^/\^yj2n. 

The results of the numerical photon-counting experiments are given in Fig. 1. For seg- 
ments of anticipated maxima and minima of intensity of the same width, I^^x, ^mm, in a region 
of the detector plate where /d(r) is much larger than the background, the visibility V is given by: 
V = (/max — /min)/(/max + Imin)- Each experiment consisted of generating photons randomly on 
the detector plate according to a probability distribution given by the numerical solution to the 
diffracted field intensity. Only the region in which /d(r) > 100[/p(r) + /pd(r)] was retained in 
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7(1 (r), which was then summed in the 2; -direction, a direction of approximate symmetry, and sub- 
sequently normalized to generate the one-dimensional probability density function used in each 
numerical trial. Numerical modelling of 10,000 successive experimental trials was performed, in 
which the visibility was repeatedly calculated for each new incident photon in the trial, with a total 
of 10,000 photons per trial. This allowed us to determine how the visibility fluctuated around the 
analytical value ofVo — 47.6%, where the fluctuations in general decreased with more detected 
photons. This produced stochastic trails for each trial, in which the largest number of photons 
where the fluctuation was greater than a given value (indicated on the horizontal axis of Fig. 1) 
was taken to be the number of photons required to reach that accuracy in the visibility. 
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